Yếu tố cường độ ứng suất là gì? Các bài nghiên cứu khoa học

Giới thiệu

Yếu tố cường độ ứng suất (Stress Intensity Factor – SIF) là đại lượng đặc trưng mức độ tập trung ứng suất ở đầu vết nứt trong vật liệu đàn hồi tuyến tính. Đại lượng này đóng vai trò quan trọng trong cơ học đứt gãy, cho phép dự đoán khả năng lan truyền vết nứt và độ bền mỏi của kết cấu chịu tải trọng tĩnh hoặc lặp.

Khái niệm SIF được phát triển từ giữa thế kỷ 20, khi George R. Irwin mở rộng lý thuyết cơ học đứt gãy để mô tả ứng xử của các vật liệu giòn như thép, gốm và hợp kim dưới tác dụng của vết nứt. Đại lượng SIF đã trở thành tiêu chuẩn vàng để so sánh độ bền đứt gãy giữa các vật liệu khác nhau và đánh giá an toàn cho công trình như cầu, bồn áp lực, và cánh tuabin.

Trong kỹ thuật, SIF giúp thiết kế kết cấu với tuổi thọ mỏi dài, giảm nguy cơ vỡ đột ngột do mỏi và ăn mòn. Việc xác định SIF chính xác hỗ trợ tối ưu kích thước vết nứt chấp nhận được, lập kế hoạch bảo trì và sửa chữa định kỳ dựa trên mô hình tăng trưởng vết nứt theo Paris Law.

Định nghĩa yếu tố cường độ ứng suất

Yếu tố cường độ ứng suất K là hằng số tỷ lệ giữa biên độ trường ứng suất σ_ij tại gần mũi vết nứt và căn bậc hai bán kính r quanh mũi nứt. Trong tọa độ cực (r,θ) tâm tại đầu vết nứt, thành phần σ_θθ có dạng:

$\sigma_{\theta\theta}(r,\theta) \approx \frac{K}{\sqrt{2\pi r}} f(\theta)$

trong đó f(θ) là hàm góc phụ thuộc chế độ tải. Giá trị K phụ thuộc hình học vết nứt, điều kiện biên và chế độ tải (mở, trượt hoặc xé).

Ký hiệu K_I, K_II và K_III lần lượt ứng với:

• K_I: chế độ I – mở vết nứt (tách hai bờ nứt).
• K_II: chế độ II – trượt trong mặt phẳng nứt (bờ nứt trượt ngang).
• K_III: chế độ III – xé vết nứt (trượt vuông góc mặt nứt).

Ba chế độ dịch chuyển vết nứt

Mode I (Opening mode) là chế độ phổ biến nhất trong phân tích vết nứt, khi lực kéo mở hai bờ nứt. Ứng suất chủ yếu sinh ra song song với mặt phẳng kéo và tập trung mạnh tại mũi nứt, dễ dẫn đến lan truyền vết nứt theo hướng vuông góc.

Mode II (Sliding mode) xảy ra khi ứng suất cắt kích hoạt trượt hai bờ nứt trong mặt phẳng bản thân vết nứt. Tải trọng cắt tạo ra moment xoắn quanh trục vết nứt và làm cho bờ nứt trượt ngang, thường gặp ở kết cấu chịu tải trọng lệch tâm.

Mode III (Tearing mode) phát sinh khi lực cắt tác động vuông góc với mặt nứt, tạo sự xé ngang hai bờ. Hiện tượng này ít phổ biến hơn nhưng quan trọng trong phân tích mỏi đứt ở vật liệu dạng tấm mỏng hoặc cấu kiện chịu xoắn nặng.

Chế độ	Ký hiệu	Chuyển vị đặc trưng	Ứng suất tập trung
Mở	KI	Δuy	σyy
Trượt	KII	Δux	σxy
Xé	KIII	Δuz	σxz

Công thức và tính toán

Đối với mẫu tấm mỏng chứa vết nứt trung tâm dài 2a, dưới ứng suất đơn hướng σ, hệ số hình học Y xác định qua hàm xấp xỉ sau:

$Y = 1.12 - 0.23\frac{a}{W} + 10.55\left(\frac{a}{W}\right)^2 - 21.72\left(\frac{a}{W}\right)^3 + 30.39\left(\frac{a}{W}\right)^4$

Do đó hệ số cường độ ứng suất chế độ I được tính:

$K_I = \sigma \sqrt{\pi a}\,Y$

với W là chiều rộng tấm. Công thức này đảm bảo độ chính xác cao khi tỉ lệ a/W nằm trong khoảng 0.2–0.8 theo chuẩn ASTM E399.

a/W	Y tương ứng
0.2	1.12
0.4	1.03
0.6	0.89
0.8	0.65

Các phần mềm phần tử hữu hạn (FEA) cũng được sử dụng để tính K thông qua J-integral hoặc Virtual Crack Closure Technique, cho phép phân tích chính xác cả chế độ pha trộn và hình học phức tạp.

Phương pháp xác định

Xác định hệ số cường độ ứng suất (SIF) có thể thực hiện theo các phương pháp thực nghiệm và số học:

• Thí nghiệm mẫu CT (Compact Tension): Mẫu CT với vết nứt chuẩn được tải trong máy kéo, theo dõi lực F và mở vết nứt Δa, sử dụng công thức:
$K_I = \frac{F}{B \sqrt{W}}\,f\!\bigl(\tfrac{a}{W}\bigr)$
• Thí nghiệm ba điểm uốn (Three-point Bend): Mẫu SENB (Single Edge Notch Bend) đo lực tải P và biến dạng khe nứt, áp dụng chuẩn ASTM E1820 để tính K_I dựa trên hệ số hình học Y.
• Phương pháp J-integral: Tính tích phân J trên đường bao quanh mũi nứt, liên hệ với K qua:
$J = \frac{K_I^2}{E'}\quad \bigl(E' = \tfrac{E}{1-\nu^2}\bigr)$
• Virtual Crack Closure Technique (VCCT): Ứng dụng trong phần tử hữu hạn (FEA), tính năng lượng cần thiết để mở rộng khe nứt nhỏ Δa và suy ra K_I.
• Photoelasticity: Mô hình trong vật liệu quang mỏng (epoxy), quan sát vân ứng suất dưới ánh sáng phân cực, xác định K từ độ giãn vân quanh mũi nứt.

Yếu tố ảnh hưởng

Giá trị SIF chịu tác động bởi nhiều yếu tố, cần cân nhắc khi tính toán và thiết kế:

• Hình học vết nứt: độ dài a, độ dày B, tỉ lệ a/W; điều chỉnh hệ số hình học Y.
• Hình dạng kết cấu: lỗ khoan giảm ứng suất (notch), góc vát mũi nứt, vết nứt cong hoặc xuyên tấm.
• Điều kiện biên: tải trọng tĩnh, tải trọng mỏi, tải trọng nhiệt độ và áp suất nội/ngoại.
• Vật liệu: mô đun đàn hồi E, hệ số Poisson ν, giới hạn đàn hồi σ_y và tính dị hướng (composite, hàn).
• Ứng suất dư (Residual stress): do hàn, gia công; cộng vào ứng suất ngoài, làm tăng hoặc giảm K_I.
• Ăn mòn và môi trường: hydrogen embrittlement, stress corrosion cracking (SCC) tăng tốc độ lan truyền vết nứt và giảm K_IC.

Yếu tố	Ảnh hưởng chính	Ghi chú
Độ dài vết nứt a	K tăng ∝√a	Quan trọng nhất trong LEFM
Độ dày tấm B	K giảm khi B tăng	Chế độ plane strain ưu tiên
Ứng suất dư	Ktotal=Kapplied±Kresidual	Phải đo hoặc mô phỏng
Ăn mòn	Tăng tốc mỏi, giảm KIC	Cần thử nghiệm SCC

Ứng dụng trong kỹ thuật

Sử dụng SIF trong thiết kế và bảo trì kết cấu công nghiệp, giao thông và năng lượng:

1. Đánh giá tuổi thọ mỏi: Tăng trưởng vết nứt theo Paris Law:
$\frac{da}{dN} = C\,(\Delta K)^m$

trong đó C, m là hằng số vật liệu, ΔK biên độ SIF lặp.

2. Khalt được khuyết (Notch design): Tối ưu hóa góc vát mũi nứt, lỗ khoan giảm ứng suất để tăng K_IC hiệu dụng.
3. Bảo trì dự đoán (Predictive Maintenance): Giám sát vết nứt qua siêu âm hoặc bề mặt (AE – Acoustic Emission), xác định ΔK và dự báo thời điểm cần sửa chữa.
4. Thiết kế vật liệu tiên tiến: Composite sợi, vật liệu gradient để giảm độ tập trung ứng suất và nâng K_IC.

Giới hạn và thách thức

Mô hình LEFM và SIF chỉ áp dụng khi vùng nhựa quanh mũi nứt nhỏ so với a (r_p≪a). Vật liệu dẻo cao, polyme và kim loại gia công bằng nhiệt không phù hợp, cần sử dụng Elastic–Plastic Fracture Mechanics (EPFM) với J-integral hoặc CTOD.

• Mixed-mode fracture: Kết hợp KI, KII, KIII khiến phân tích phức tạp; cần mô hình số học đa chế độ.
• Ảnh hưởng microstructure: Hạt, vết nhựa, kẽ hở ảnh hưởng cục bộ, gây dao động K và không đồng nhất.
• Độ chính xác thí nghiệm: Đo a và lực F nhỏ, sai số lốc kéo đòi hỏi thiết bị chuẩn và kỹ thuật viên giàu kinh nghiệm.

Hướng nghiên cứu tương lai

Ứng dụng công nghệ không dây và IoT để đo real-time SIF qua cảm biến strain gauges in situ, truyền dữ liệu về trung tâm theo dõi (NIST).

Sử dụng machine learning và công nghệ big data phân tích chuỗi dấu hiệu AE, siêu âm và tải trọng, dự báo tăng trưởng vết nứt với mô hình tiên tiến (ScienceDirect).

Phát triển vật liệu tự liền vết nứt (self-healing composites) có khả năng giảm hiệu quả ΔK thông qua mạng lưới polymer và vi nang chứa chất chữa lành.

Danh mục tài liệu tham khảo

• Anderson T.L. “Fracture Mechanics: Fundamentals and Applications.” CRC Press, 2017.
• Paris P., Erdogan F. “A Critical Analysis of Crack Propagation Laws.” J. Basic Eng. 1963;85(4):528–534.
• ASTM E399-20. “Standard Test Method for Linear-Elastic Plane-Strain Fracture Toughness of Metallic Materials.” ASTM International.
• Elber W. “Fatigue Crack Closure Under Cyclic Tension.” Eng. Fract. Mech. 1970;2(1):37–45.
• NIST. “Fracture Mechanics Handbook.” nist.gov.
• Rice J.R. “A path independent integral and the approximate analysis of strain concentration by notches and cracks.” J. Appl. Mech. 1968;35:379–386.